Polar and non-polar atomic motions in the relaxor ferroelectric PLZT from dielectric, 

anelastic and NMR relaxation 
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Dielectric, anelastic and La NMR relaxation mea- 
surements have been made on the relaxor ferroelectric 
Pbi-3;r/2La:,Zro.2Tio.803 (PLZT) with x = 0.22. The dielec- 
tric susceptibility exihibits the frequency dispersive maximum 
due to the freezing of the polar degrees of freedom around 
T ~ 250 K. The anelastic and especially NMR relaxation, 
besides this maximum, indicate an intense and broad compo- 
nent at lower temperatures, attributed to rotational modes of 
the O octahedra, weakly coupled to the polar modes. It is dis- 
cussed why such short range rotational instabilities, known to 
occur in the Zr-rich rhombohedral region of the PLZT phase 
diagram, might appear also in the Ti-rich region. 



I. INTRODUCTION 

Much effort is being spent during the last years in 
trying to understand the so-called relaxor state of dis- 
ordered ferroelectrics, and notably the perovskite relax- 
ors, with general formula ABO3. The technological in- 
terest arises from the fact that such materials may be 
prepared with very high values of dielectric, piezoelec- 
tric, electrostrictive and electro-optical constants in wide 
temperature ranges, finding applications like non- volatile 
and DRAM memories, capacitors, microelectromechan- 
ical devices, electro-optic modulators, optical switches 
etc. 

Such materials have disorder in the cation sublattices 
that inhibits the formation of long range ferroelectric 
order, ^ and present a phenomenology similar to that 
found in spin glasses^ and other glassy systems, with 
freezing of the polar fluctuations, separation of field 
cooled from zero field cooled susceptibility curves'^ and 
non equilibrium phenomena like aging and memory.'''^ 
After the early explanation of the broadened maximum 
of the dielectric susceptibility in terms of a distribu- 
tion of Curie temperatures, according to the local com- 
position, other models have been proposed, e.g. in 
terms of quenched random fields due to compositional 
fluctuations,^ and spin-glass state induced by interac- 
tions among the polar clusters^ or random fields.^ Re- 
cently, the spherical random-bond random-field model* 
has been developed, as a useful framework to understand 
several properties of the relaxor ferroelectrics. 

Regarding the microscopic mechanism of the elec- 
tric polarization, the focus has been on perovskites 



having charge disorder in the B sublattice, like 
Pb(Mgi/3Nb2/3)03 (PMN). In these systems detailed 
studies of the atomic displacements have been carried out 
by neutron^'^" and synchrotron x-ray-'^^ diffraction, and 
by x-ray absorption fine structure, while the relevant 
phonon instabilities have been studied by inelastic neu- 
tron spectroscopy. ^"^"^^ Much recent experimental^^'-'^^ 
and theoretical^* work is devoted even to the well 
known ferroelectric system PbZri_yTij,03 (PZT) near 
the boundary at y ~ 0.48 between the Zr-rich rhombo- 
hedral and Ti-rich tetragonal phases. 

Another perovskite relaxor system used in several ap- 
plications is (Pb/La)(Zr/Ti)03 (PLZT), obtained from 
PZT by substituting Pb^+ with La'^^. Since in this case 
the charge disorder is in the A sublattice, many of the 
considerations valid for B-disordered relaxors do not hold 
for PLZT. The most explored region of the PLZT phase 
diagram is the one around 65% of Zr, while only few 
studies are devoted to the Ti-rich one,"'^^"^^ which show 
particularly high quadratic electro-optic effect and elec- 
trostrictive strain, useful for quadratic electro-optic mod- 
ulators, optical phase retarders, electro-optical shutters 
and electrostrictive actuators. 

In this report we present the results of a study of the 
atomic dynamics accompanying the freezing of polariza- 
tion in La- and Ti-rich PLZT by measuring its dielectric 
susceptibility, elastic compliance and NMR relaxation 
rate of the ^^^La nuclei. These experiments probe the 
fluctuations of the electric polarization, of strain and of 
the environment of the La nuclei, respectively. 



II. EXPERIMENTAL 

Ceramic (Pbi_3a;/2La^)(Zro.2Tio.8)03 
(PLZT 100a::/20/80), with the charge compensating va- 
cancies in the Pb sublattice, has been prepared by using 
the mixed-oxide method according to the following pro- 
cessing route. The oxide powders were wet ball milled 
with zirconia balls in the stoichiometric amount for 24 h, 
the suspension was then freeze dried, sieved to 200 ^m 
and calcined at 850 °C for 4 hours. Samples of the cal- 
cined powders with x = 0, 2, 12 , 22% were pressed into 
bars and sintered at 1250 °C for 2 h, packed with PbZr03 
-|- 5wt% excess Zr02 in order to maintain a constant PbO 
activity at the sintering temperature. 
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Phases and microstructurc were investigated by x-ray 
diffraction (XRD) analysis and by SEM on fracture on 
a series of samples with < a; < 0.22. For x = 0.22, 
the average grain size determined from SEM was 3.2 /xm, 
and the density was about 95% of the theoretical one. 
The XRD patterns in Fig. 1 show the presence of pure 
perovskitic phase with symmetry varying from tetrago- 
nal to pure cubic as the addition of La increases; the 
tetragonal distortion calculated over appropriate peaks 
of the diffractograms is shown in Fig. 2 to vary almost 
linearly from 5% for PZT 20/80 to < 0.2% for PLZT 
22/20/80. The composition under investigation here is 
PLZT 22/20/80, with a cubic cell and without long range 
ferroelectric order; in fact, ferroelectric hysteresis mea- 
surements showed a very slim loop with a remanent po- 
larization of about 1 iiC/cm^ (typical values in the ferro- 
electric state are 30 /zC/cm^). Samples were cut as thin 
bars approximately 45 x 4 x 0.5 mm"^, the shape required 
by the anelastic spectroscopy experiments. After cutting, 
the samples where annealed for 15 h in air at 750 °C; the 
electrodes for the anelastic and dielectric spectroscopy 
measurements were applied with silver paint. 

The dielectric susceptibility x = x' + *x" was measured 
with a HP 4194 A impedance bridge with a four wire 
probe and a signal level of 0.5 V/mm, between 200 Hz 
and 1 MHz. The measurements were made on cooling at 
1-1.5 K/min from 570 to 200 K. 

The complex Young's modulus E (w, T) = E'+iE" was 
measured on cooling by suspending the bar on thin ther- 
mocouple wires and electrostatically exciting its 1st and 
5th flexural modes, whose resonance frequencies where 
Wi/2TT ~ 1.5 and 20 kHz. The real part of the Young's 
modulus i s relat ed to the resonance frequency through 
uji = OLi^ E' / p, where is a geometrical factor and 
p the mass density. The elastic energy loss coefficient, 
or the reciprocal of the mechanical quality factor, is 
g-i {uj,T) = E" IE' = a" /s', where s ^ s' + is" = E-^ 
is the elastic compliance; the Q^^ was measured from the 
decay of the free oscillations or from the width of the res- 
onance peak. The elastic compliance s is the mechanical 
analogue of the dielectric susceptibility x- 

The ^'^'^La NMR spectra and spin-lattice relaxation 
rates 2W have been obtained by standard pulse tech- 
niques, operating in the quadrupole perturbed Zeeman 
regime, with external magnetic field Hq varying from 1.8 
to 9 Tesla. The temperature range explored was 120—600 
K. The spectra were obtained by the Fourier transform 
(FT) of the half of the echo signal and also by the enve- 
lope of the echo under sweeping the RF frequency. The 
echo signal arises from the powder distribution of the 
second-order quadrupole shift of the central line. From 
FT the central component showed a width of about 65 
kHz. The satellite components are spread over a wide 
frequency range because of the first-order perturbation 
effects due to the electric quadrupole interaction. The 
pulse length maximizing the ^'^^La quadrupolar echo sig- 
nal was found about a factor 2 smaller than the pulse 
expected for La nuclei in the absence of quadrupole in- 



teraction, indicating that the satellite lines are almost 
not irradiated. The relaxation time measurements were 
performed by monitoring the recovery of the echo ampli- 
tude after irradiation of the central line with a sequence 
of saturating pulses in a time much smaller than the re- 
covery time. Even though the acquisition of the recovery 
of the nuclear magnetization lasted several hours, aging 
did not affect the measurements, since the aging effects 
decrease with increasing frequency. In order to assess the 
lack of influence from aging, the recovery of the nuclear 
magnetization was recorded using different sequences of 
delays and checking that the amplitude of the signals did 
not depend on the chosen sequence. 

For the sake of illustration of the relationship between 
the different experimental techniques used in the present 
work, we anticipate here that in the ideal case that the 
main contribution to the nuclear relaxation process orig- 
inates from the same type of atomic displacements u {t) 
that determine the dielectric and elastic uniform suscep- 
tibilities, 2W could be qualitatively written in the form^^ 

2W = A^ j {u (t) u (0)) e-'^'dt = £■ J (w) (1) 

where ^ is a constant including the quadrupole mo- 
ment of the ^■^^La nucleus and the derivatives of the 
electric field gradient (EFG) components at the nuclear 
site with respect to the displacement 7i, while J (uj) is 
the spectral density of the correlation function for u{t). 
In this model situation, a simple relationship between 
the NMR relaxation rates and the dielectric data would 
be possible through the fluctuation-dissipation theorem, 
x" (cj) cx ^J(cj). Later on we shall discuss the break- 
down of such a direct connection. 



III. RESULTS 

Figure 3 presents the real and imaginary parts 
of x('^)^) S'lid s(w,T), together with one of the 
Xnmr('^'^) curves defined as Xnmr('^'^) = 
according to the aforementioned simplified model (for 
w = 27r X 54 MHz). Even though all three types of re- 
laxations present a maximum in correspondence with the 
relaxor transition around 250 K, there are substantial dif- 
ferences in the position and frequency dispersion of this 
peak, and in its low temperature tail. The comparison 
between dielectric susceptibility and elastic compliance 
is more direct, since both involve uniform susceptibili- 
ties X (? = 0, Lo) in the same frequency ranges. The dif- 
ferences between the x and s curves indicate the pres- 
ence of additional relaxation modes, besides the freez- 
ing of the polar degrees of freedom probed by x- For 
Xnmr ("-^i ^) on® should take into account a sizeable con- 
tributions from dispersive relaxation modes at different 
wave vector which are not included in the uniform sus- 
ceptibilities. 
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A. Dielectric susceptibility 



In T2 = In TO + E/ [T - Tvf) 



(5) 



The x' {^^ T) curves in Fig. 3 are in good agreement 
with the series of x' (t^, T) curves for PLZT a;/20/80 with 
16 < X < 20 aheady reported, the maximum being 
sHghtly smaller and about ^ 50 K lower in temperature 
with respect to the highest x reported there. In order 
to make a comparison with the NMR 2T4^ curves, we de- 
termine a fitting expression for x(w,r) that can be ex- 
trapolated at the NMR frequencies. The interpolation of 
the dielectric susceptibiUty involves the choice of a static 
susccptibihty x (0, T), and of a frequency dispersion fimc- 
tion, generally taken as a superposition of elementary re- 
laxation functions with a suitable distribution g (In r) of 
relaxation times, so that the dielectric susceptibility is 
expressed as 



with To = IQ-^" s, E = 1450 K and Tyr = 190 K, 
while the fastest relaxation time has been set to ti = 
5 X 10~^^ s. The x (<^, T) curves have been obtaining by 
numerically integrating cq. (2) and arc shown in Fig. 5. 
The fit is optimized on x" (w, T) at the highest frequency, 
in order to obtain a reliable extrapolation of the spectral 
density 



J(w,T)cx -x" (w,T) 



(6) 

at the NMR frequencies. For w = 0, x" 
has the analytical form x" = x{^-,T) (\nT2/Ti)~^ 
[arctan (ti)T2) — arctan (wti)], which is valid also for w ^ 
and T below the maximum, where the shape of the cut- 
off at ijJT2 '> 1 has no influence. The low T approxima- 



X (w, T) = x{0,T) J dluT g (In r) • (2) tion is therefore x" (w, T) ~ x (0, T) f (^^) , where 



In relaxor ferroelectrics, x (0, T) follows a modified Curie- 
Weiss law at temperatures sufficiently far from the re- 
laxor peak, but there is no divergence on approaching it. 
For the present purposes, we limit ourselves to finding a 
phenomenological expression for x (0, T) that is suitable 
to fit the data down to temperatures as close as possible 
to the relaxor peak. An excellent fit can be obtained with 
a logistic function, which saturates at low T instead of 
diverging. 



x(0,T) = xo + 



Ax 



1 + {T/Tof 



(3) 



with xo = 270, Ax = 9260, Tq = 319 K, p = 5.65. The 
value of Xo seems too high for representing the infinite 
frequency limit Xoc, and is therefore considered as a fit- 
ting parameter in (3), while Xoo is neglected; the fitting 
curve is shown in Fig. 4. 

Regarding the frequency dispersion, we tried with 
exponential^^ and uniform^^"^^ distribution functions 
in g{\nT). The first type of distribution function 
was adopted to interpret dielectric measurements on 
Pb(Mgi/3Nb2/3)03,^'' but could hardly reproduce the 
present data. Instead, we could interpolate reasonably 
well both x' and x" with a modified uniform distribution 
between Inn and lnr2 (T), with the cutoff at the longest 
time broadened over w (T) in the logarithmic scale: 



gilnr) 



goi[l-tanh(ilii^)] forr>Ti 



for T < Ti 



(4) 



50 = 2/ 



lnT2/Ti -|- win 



2coshfi^i^ 



Tl 



with w (T) = 0.15 X (lnT2 (T) - Inn) scaling with the 
distribution width. The longest relaxation time at half 
maximum of the distribution has been assumed to fol- 
low the Vogel-Fulcher law, similarly to what observed for 
other compositions of PLZT^'^^ 



E2-E1 y' 

the linear increase with T above Tvf comes from the T 
dependent normalization factor (lnT2/n) therefore, 
considering the weak temperature dependence of x (0, T) 
below the maximum, Tvf is the temperature at which the 
linear extrapolation of x" below the maximum vanishes. 
Of course, there is no real divergence of the distribution 
width and the dielectric susceptibility does not vanish 
below Tvf, but the Vogel-Fulcher hypothesis allows the 
dispersive maximum to be well described. 



B. Elastic compliance 

The elastic compliance curves s (w, T) measured at the 

1st and 5th resonance frequencies (1.5 and 20 kHz) are 
shown in Fig. 3 together with the susceptibility curves. 
The curves are normalized to the extrapolated value of s' 
at infinite temperature, Sqo, which should be the equiva- 
lent of Xoo- An important difference between the dielec- 
tric and the clastic susceptibilities, however, is that in 
relaxor ferroelectrics Xoo brings a negligible contribution 
to the dielectric response; in fact, Xoo is expected^* to 
be < 10 (even though Xoo ^ 200 has ben reported in 
PLZT 9/65/3529) and here it is x' ~ 7000 at the max- 
imum. Instead, the background lattice contribution to 
s' is much more important: the peak value of s' is only 
1.6 times larger than Sqo. Since the lattice contribution 
to s' is also temperature dependent, it is easier to com- 
pare only the imaginary parts of the different suscepti- 
bilities, to which Soo and Xoo do not contribute. The 
elastic s (w, T) differs from x (w, T) in three ways: i) the 
maxima of the real and imaginary parts are at lower tem- 
perature; ii) the frequency dispersion is smaller both in 
temperature shift and change of intensity, and the am- 
plitude dependence of s' is reversed with respect to x'i 
in) the low temperature tail is much higher than in the 
dielectric case (s" (150 K) /s" (Tmax) = 0.5 for the clas- 
tic case at 20 kHz, while x" (150 K) /x" (T^ax) = 0.07 at 
10 kHz). The comparison between x" and s" is rather 
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straightforward, since they both represent uniform sus- 
ccptibihties. The only differences are in the relaxation 
strengths, which are sensitive to changes in electric dipo- 
lar degrees of freedom in one case and elastic quadrupolar 
in the other. In the simplest case of independent relaxing 
units, e.g. off-centre atoms hopping between positions i 
and j with a correlation time r one has^^ 



Ax(0,T) = 



{Apf 
knT 



As (0,T) = 



{AXf 



(7) 



where Ap = p*-*^ — p^^^ is the change of the electric dipole 
during the jump, and AA — A*^*' — A'--'-' is the change of 
the elastic quadrupolc (generally called elastic dipole^^ 
for analogy with the electric and magnetic cases). The 
relaxation of a same unit may produce a change in p 
but not in A, e.g. the reorientation of an off-centre unit 
by 180° changes p to — p but does not affect A, which 
is a centrosymmetric strain (representable as a strain el- 
lipsoid). Conversely, there may be relaxation modes af- 
fecting A but not p, e.g. rotations or distortions of the 
octahedra. The differences in the dielectric and elastic 
susceptibilities indicate that the latter is little affected 
by the dipolar relaxation, since the peak in s" is not 
proportional to the relaxor peak in in addition, there 
are other relaxation modes almost invisible to x-> but con- 
tributing to s below the relaxor transition. These modes 
are faster and only loosely follow the dipolar freezing, as 
indicated by the lower temperature of the peak and the 
huge low temperature tail. 

The role of clcctrostrictive strain in driving the 
anelastic relaxation in relaxor ferroelectrics has been 
emphasized, ^^'^"^ but cannot account for the present re- 
sults; in fact, if the strain fluctuations were simply driven 
by the dipolar ones through electrostriction, s" would 
follow x" much more closely, since the electrostrictive 
coupling is expected to be weakly dependent on temper- 
ature. 



C. NMR 

Under the conditions of measurements outlined in Sec- 
tion II, by solving the master equations for I = 7/2, the 
recovery of the echo signal is expected to be described by 
the law 



y{t) = 0.012 e 



0.714e-i-33W'* + 



(8) 



+0.068 e-2-4^*+ 0.206 e-3-8VVt 



In deriving Eq. (8) the relaxation transition probability 
Wi and W2 (Am = 1 and Am = 2 transitions respec- 
tively) have been assiimcd practically equal, as usual in 
powdered samples. -^"^ At high temperatures, well above 
the ferroelectric relaxor transition temperature, the re- 
covery of the echo signal is rather well fitted by Eq. (8), 
as shown by the solid line in Fig. 6. One notes that, 
according to Eq. (8), the relaxation rate 2W can also be 



obtained from the time t* where y{t*) = 0.41, yielding 

2W = {t*)-\ 

In the low temperature range Eq. (8) does no longer 
describe the recovery of the echo (see Fig. 6). This is 
the obvious consequence of the fact that, on approach- 
ing the relaxor transition temperature, each exponential 
in Eq. (8) becomes stretched. A meanigful fitting ac- 
cording to Eq. (8) is hardly possible, and thus in the low 
temperature range (T < 300 K) an effective 2W has been 
extracted as the inverse of the time t* . 

The ^^^La relaxation rates 21^ are reported in Fig. 
7 for three representative measuring frequencies aj/27r. 
From 700 K, 2W increases on decreasing temperature 
passing through a maximum and then decreasing again. 
The maxima in 2W occur in the temperature range where 
the maxima in imaginary part of the dielectric constant 
arc also detected (see Fig. 3). This suggests that, at 
least in the high temperature range, the mechanisms re- 
sponsible of the reorientations of the electric dipole mo- 
ments arc also involved in the ^'^^La relaxation process. 
It should be remarked that, by varying the magnetic field 
from ifo = 1-8 T to ilo = 9 T, the maxima of 21^ do 
not sizably shift in temperature. This direct qualitative 
observation by itself suggests that the freezing process 
involves a large distribution of correlation times, as can 
be expected for a strongly disordered relaxor. 

In the weak collision approach, for overdamped phonon 
modes and/or localized modes in disordered systems (di- 
rect relaxation process) the relaxation rate driven by the 
time dependence of the EFG components can be written 
in the form anticipated in Eq. (1), where the spectral 
density Ji{ui) and J2(2w) have been assumed practically 
coincident. We can explicitate the constant A"^ and write 
2W in the form 2W ~ 5(wq) fiJ{Lu), (luq) fi being an av- 
erage effective strength of the time dependent quadrupole 
coupling and J{uj) the spectral density of the reduced cor- 
relation function.^'^ For monodispersive processes, with 
a single correlation time, the spectral density would be 
J(a;) = 2r/(l + uj^t^). In the fast motions regime, i.e. 

u), one would have a frequency independent 21^ 
of the order of 2W ~ 12(wQ)//r, while for slow motions 
(t"1 < uj) it is 2W ~ 10(w^)/;(w2t)"^ From the value 
of 2W at the maximum an order of magnitude of the 
fluctuating quadrupole coupling constant can be derived: 
(a;|)}f ~ (2iy„/w)V2 ~ 50 - 100 kHz, a reasonable 
value for the motions around the equilibrium position 
of the ions. In the present case, the lo dependence of 
2W is reduced by the broad distribution of correlation 
times typical of relaxors. We first tried with the same 
uniform distribution function adopted for the dielectric 
case, but it is apparent from Fig. 8 that the extrapola- 
tion of the {T/ui) x" (w, T) curves to the NMR frequencies 
does not properly reproduce 2W. We therefore tried with 
the distribution proposed by Lu and Calvarin^^ (LC) for 
their dielectric relaxation data in Pb(Mgi/3Nb2/3)03 and 
Pb2KTa50i5. The LC distribution is assumed in the 
form 
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P{E) = 



^-l^{Ec-E)/w ioj.E>Ec 

for S < Ec 



(9) 



where E = KV is the activation energy for thermal fluc- 
tuations of a polar cluster in the paraelectric matrix, pro- 
portional to its volume. Ec = KVc is the minimum ac- 
tivation energy corresponding to a critical volume of the 
order of the volume of a crystalline cell. The temper- 
ature dependence of the distribution width w has been 
assumed of the form w = Eqc^^^'^. The correlation times 
are related to the activation energies by r = roe^/"^. In 
Fig. 7 the lines represent the temperature dependence of 
2W numerically obtained in the assumption of LC dis- 
tribution with To = 10~i3 = iQQ Eo = 115 
K and Ei = 750 K. The flt of the experimental results, 
although better than the simple extrapolation of the di- 
electric ones, is again poor. 

A detailed scries of measurements as a function of lu 
has been carried out at constant temperature, T = 140 
K, below the maxima in 2W vs T. In Fig. 9 the experi- 
mental data at T = 140 K have been compared with the 
frequency dependence expected on the basis of the re- 
laxation rates according to the LC distribution (dashed 
line) and to a uniform distribution (dotted line). It is 
evident that both the absolute values as well as the fre- 
quency dependence do not agree with the theoretical ex- 
pectations according to the above models. Here we note 
that a good fit could be obtained if a frequency indepen- 
dent contribution, arising from a second type of motion 
efifective in the relaxation process, is added to the fre- 
quency behavior of the LC distribution (solid line). As 
noted above, an almost frequency independent contribu- 
tion implies a prevalence of fast motions, with lot <C 1 
still at 140 K, consistent with the broad contribution at 
low temperature both in 2W and in the anelastic s". We 
conclude that the ^''"La NMR relaxation contains addi- 
tional contributions with respect to the dipolar freezing 
dominating the dielectric susceptibility, especially in the 
low temperature region. 



IV. DISCUSSION 

A. Non-polcir relaxation modes 

The comparison between dielectric, anelastic and NMR 
relaxations indicates the existence of at least two distinct 
types of atomic motions: i) polar modes that dominate 
the behavior of giving rise to the typical relaxor max- 
imum, but are less evident in s" and 2W: ii) other modes, 
affecting much less x and dominating the nuclear mag- 
netic and anelastic relaxations. The second type of modes 
prevails at lower temperature and should be essentially 
non-polar, since it does not appear in the dielectric re- 
laxation. 

A comparison between dielectric and anelastic relax- 
ation in PLZT 9/65/35, shows even larger differences,^^ 
since s" completely lacks the relaxor peak and only 



presents a step-like increase at a temperature ^ 30 K 
lower than that of the relaxor peak in x"; this is shown 
in Fig. 10. It was proposed that the main difference be- 
tween the two types of responses consists in the lack of 
sensitivity of the elastic compliance to the 180° switching 
of the polarization, which would be the main responsible 
for the relaxor peak; the non-180° polarization dynamics, 
contributing also to s", would instead have an increasing 
role at lower temperature.''^ The NMR data, however, 
suggest a different picture; in fact, there is no reason 
to suppose that the ^^®La nuclear relaxation should be 
much more sensitive to non-180° rather than to 180° po- 
larization dynamics, and therefore the explanation of the 
predominant weight of relaxation in 2W at temperature 
lower than that of the relaxor peak can hardly rely on 
non-180° polarization dynamics. Rather, this indicates 
that the NMR relaxation probes a non-polar dynamics, 
not appearing in presumably the same type of mo- 
tion appearing in the anelastic response. Motions of this 
type may be rotations or distortions of the octahedra and 
antiferrodistortive displacements of the cations. 



B. Rotations of tlie octaliedra 

Rotations or deformations of the octahedra would cer- 
tainly provide a strong contribution to 2W , since a La ion 
is surrounded by 12 nearest neighbor O atoms at a/\/2, 
and 8 second neighbor Ti/Zr atoms at \/3/2a; the con- 
tribution of the atomic shifts to the EFG fluctuations 
decreases with the 6th power of the distance from the 
La site, so that the contribution from motions of the O 
octahedra (or of La with respect to the surrounding octa- 
hedra) is enhanced by 5 times over the contribution from 
the motion of the Zr/Ti atoms. Since the microscopic po- 
larization in PLZT is mainly attributed to the off-centre 
displacements of the cations with respect to the O octa- 
hedra, and there is little correlation between these polar 
fluctuations and 2W, it seems reasonable to assign the 
main contribution to the NMR relaxation to rotations or 
deformations of the octahedra. 

Distortions of the octahedra of the Jahn- Teller type 
are found in perovskites with cations supporting mixed 
valence, like manganites, but this is not the case of ferro- 
electric PLZT, where both Zr and Ti have fixed valence 
+A. Jahn- Teller polarons in ferroelectric perovskites may 
be observed only after suitable doping, e.g. in Nb- 
doped BaTiOa where each Nb'''+ is compensated by one 
Jahn- Teller active Ti"'+; our samples might well contain 
some minority charged defects, like O vacancies, but they 
are insulators and a mechanism like a Jahn- Teller polaron 
cannot be considered as source of the intense relaxation 
observed below the relaxor transition. 

Rotations of the octahedra are instead typical of all 
perovskites, including the ferroelectrc ones, and the low 
temperature phases of Zr-rich PZT are the result of the 
condensation of unstable rotational modes of the O octa- 
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hedra, together with fcrro- and antifcrroclcctric modes of 
the cation sublatticcs.^^'^'' The PZT phase diagram^^^'^^ 
grossly consists of a cubic high temperature phase, which 
transforms into ferroelectric tetragonal for a Ti fraction 
y > 0.48 and ferroelectric rhombohedral below the al- 
most temperature independent boundary y ~ 0.48. In 
the tetragonal phase, the A (Pb) and B (Zr/Ti) sub- 
lattices are shifted along the c axis with respect to the 
O octahedra, which are additionally elongated. This is 
seen as the result of the condensation of the unstable 
Fis phonon. In the rhombohedral phase the cations are 
shifted along the pseudocubic [111] direction, and, be- 
low a certain temperature, the octahedra rotate about 
the polar [111] axis, giving rise to the low-temperature 
rhombohedral phase. Additional rotational instabilities 
of local character have been proposed to exist in Zr- 
rich PZT^s and especially in PLZT.35'37 viehland and 
coworkers found | (HI) and ^ (HO) superlattice spots 
in electron diffraction by TEM, which they associate 
with coordinated rotations of the octahedra by an angle 
±(j) about the pseudocubic c axis, with successive lay- 
ers in phase (condensation of the M3 phonon mode^^ 
producing ^ (HO) reflections) or antiphase (R25 mode 
producing 5 (111) reflections). These reflections progres- 
sively disappear with increasing Ti content, and in PLZT 
12/40/60 only weak spots were found in a minority frac- 
tion of the selected areas. "^^ It has later been argued^® 
that such superlattice reflections might mainly arise from 
additional shifts of the Pb atoms along pseudocubic [100] 
directions, coupled with the octahedral rotations; in ad- 
dition, since these reflections are not observed by x-ray 
or neutron diffraction, the regions of additional M3 or 
i?25 ordered rotations might exist only near the surface 
of grains, at least in La-free PZT. Rotational instabilities 
of M3 and i?25 type have also been found in flrst-principle 

calculations of the phonon instabilities in rhombohedral 
P2T.34 

While the existence of these rotational instabilities is 
well documented in the rhombohedral region of the phase 
diagram of PZT and PLZT, we are not aware of any 
such observation in Ti-rich PLZT. Therefore, octahedral 
rotations seem at first unlikely to occur in our PLZT 
22/20/80, considering that the parent PZT 20/80 is well 
within the tetragonal region of the phase diagram. We ar- 
gue, however, that the substitution of Pb with La changes 
the situation considerably with respect to PZT. First of 
all, XRD measurements on a series of samples with in- 
creasing La content (Fig. 1) show that the addition of 
La suppresses the tetragonal distortion, which is 5% in 
PZT 2/20/80, 2% in PLZT 12/20/80, and < 0.2% in 
PLZT 22/20/80, implying the absence of a strong sta- 
bilization of the tetragonal phase in the present case. A 
similar behavior has been found in PLZT a;/40/60.^'' The 
La addition should also favor the rotational instability at 
the expenses of the tetragonal one. In fact, calculations 
of the phonon dispersions^^'^s gjiow that both the Fis 
tetragonal and the rotational R25 instabilities exist in 



PbTiOa and PbZrOa, although the ri5, mainly consist- 
ing in a displacement of Pb along (001),'^^ is stabilized in 
PbTiOs by the tendency of Pb to shift in that direction 
in order to hybridize and form 4 short Pb-0 bonds; this 
is possible, since also the Ti ions shift off-centre from the 
octahedra in the same direction, away from the O atoms 
with short Pb-0 bonds. In PbZrOs instead, where the 
ideal Zr-0 bonds are longer than the Ti-0 ones, the ro- 
tational instability prevails in order to accommodate the 
larger octahedra; this fact is usually expressed in terms of 
the tolerance factor*^ t = (va + vq) /\/2 (re + ro) which 
is 1 if the ideal B-0 distance matches the ideal A-0 dis- 
tance, and is < 1 when the octahedra are too large to fit 
the cell and tend to rotate to maintain their size. The 
substitution of Pb with La should favor the rotational 
instability because: i) it decreases the tolerance factor, 
due to a smaller ideal La+'^-O bond length (2.76 A) with 
respect to Ph'^+-0 (2.89 A),"*! and in fact XRD showed 
that the cell volume passed from 65 A^ in PZT 20/80 to 
63.4 A3 in PLZT 22/20/80 (Fig. 2); ti) L&^+ does not 
have the ability of stabilizing the ferroelectric tetragonal 
mode by a large electronic hybridization with O, as Pb 
does.18.41 

We think that the best candidates for the non-polar 
motions we observe are rotations of the octahedra related 
to the antiferrodistortive R25 mode described above. 
First-principle calculations of the phonon dispersions 
showed that such a mode is ^instable also in PbTiOs, 
although not enough to produce the ground state. ^® The 
dispersion curves indicate a strong correlation of the co- 
operative rotations about c within the ab plane, due to 
the corner sharing of the octahedra, but little correlation 
along c (the same mode, with full correlation also along 
c, is responsible for the cubic to tetragonal transforma- 
tion in SrTiOa"*^). Due to the lack of correlation between 
different planes even in pure PbTiOs and to the strong 
disorder in the A sublattice in PLZT 22/20/80, the lo- 
cal condensation of such a mode over only few octahedra 
should be possible, without producing superlattice peaks 
in diffraction patterns. In addition, the small correla- 
tion length and frustration from disorder would make it 
possible to switch between differently rotated configura- 
tions. Switching between rotated and non-rotated con- 
figurations would be accompanied by a change of the in- 
plane shear strain (corresponding to A A 7^ in Eq. (7)) 
and therefore by anelastic relaxation, while the NMR re- 
laxation would be sensitive also to switching between cor- 
related rotations by opposite angles (which leave strain 
unchanged) and to the anharmonic dynamics related to 
such an unstable mode. 

Finally, the rotations of the octahedra are only weakly 
coupled to the cation displacements in PLZT,*^ the ma- 
jor coupling perhaps coming from the tendency of Pb 
to form four short Pb-0 bonds. In fact, recent powder 
neutron diffraction experiments on rhombohedral PZT 
have been interpreted assuming that, besides the estab- 
lished ferroelectric displacements of Pb along the pseu- 
docubic (111) direction, additional (100) shifts of 0.2 A 
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are induced by the tendency of Pb to form the fourfold co- 
ordination with the surrounding O atoms, with the choice 
between x, y or z dictated by the local tilt pattern of the 
octahedra.^^ Then, the weak coupling between the ferro- 
electric cation displacements and the octahedra rotations 
justify the fact that the latter are still fast at tempera- 
tures at which the polar modes are completely frozen out; 
on the other hand, the weak low temperature tail found 
also in the dielectric x" might result from the supposed 
additional degree of freedom of Pb,^^ which is coupled 
to the octahedra (besides the fact that disorder might 
induce a polar component also to the rotations of the 
octahedra). 

Further support to the hypothesis that rotations of 
the octahedra are responsible for the low temperature 
NMR and anelastic relaxation comes from the compari- 
son of the anelastic spectra of PLZT 22/20/80 and PLZT 
9/65/35,^^ shown in Fig. 10. The latter composition is 
in the rhombohedral region of the phase diagram, where 
the rotational instabilities are documented, and in fact 
the low temperature component alone determines the 
anelastic spectrum of PLZT 9/65/35, whereas it is less 
important in PLZT 22/20/80. 

C. Antiparallel cation motion 

In principle, also correlated antiparallel motions of the 
A or B cations might explain the discrepancies between 
the dielectric relaxation on one side and the NMR and 
anelastic relaxation on the other side, since they produce 
fluctuations both in shear strain and on the EFG at the 
La site, but negligible dielectric response. Such modes 
are responsible for the antiferroelectric phases in the Zr- 
rich end of the PLZT phase diagram, but we do not see 
reasons why they should have importance in PLZT with 
80% of Ti. In addition, such antiferroelectric modes in- 
volve displacements of the Pb atoms;^^ therefore, their 
condensation should be hindered by the strong disorder 
in the A sublattice, and, if it occurred, it would result in a 
polar mode, due to the charge variations among different 
A sites (Pb^+j La'^"'" and vacancies). We conclude that 
the extra contribution to the relaxation process found 
in the anelastic and NMR experiments is rather due to 
modes involving the O octahedra. 

V. CONCLUSION 

The combination of dielectric, NMR and anelastic 

spectroscopics in the rclaxor ferroelectric PLZT reveals 
that, besides the polar fluctuations freezing at the re- 
laxor transition, other fluctuations of the atomic posi- 
tions exist, that are faster and weakly correlated with 
the polar ones. In fact, the imaginary dielectric suscepti- 
bility x" (w, T) presents the usual peak, with a frequency 
dispersion signaling the freezing of the polar degrees of 



freedom, and a weak tail at lower temperature, while the 
curves of the imaginary elastic compliance s" {oj,T) and 
especially of the ^^^La NMR relaxation rate 2W versus 
temperature present a more intense component at lower 
temperatures. This low temperature relaxation has been 
assigned to non polar modes, which have been argued 
to correspond to rotations of the octahedra; in partic- 
ular, on the basis of the calculated phonon dispersion 
curves in the PZT system, it is suggested that the 
condensation of the R25 mode might lead to clusters of 
rotated octahedra about the pseudocubic c axis. The 
correlation lengths for such correlated rotations should 
be very short, due to the high concentration of La'^^ and 
Pb vacancies, which hinder the otherwise predominant 
ferroelectric tetragonal distortion. The small correlation 
length and the weak coupling with the other ferroelectric 
modes would make it possible for such clusters of octahe- 
dra to switch between differently rotated conflgurations, 
giving rise to the observed low temperature relax;ation. 
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VII. FIGURE CAPTIONS 



FIG. 1. Diffractograms of PLZT x/20/80 at four La con- 
tents. 

FIG. 2. Vanishing of the tetragonal distortion with increas- 
ing La content. 

FIG. 3. Real (upper panel) and imaginary parts (lower 
panel) of the dielectric susceptibility x (right ordinates) and 
elastic compliance s (left ordinates), measured at the frequen- 
cies indicated in the figure. Also shown is Xnmr defined in 
the text measured at 54 MHz. 

FIG. 4. Real part of the dielectric susceptibility measured 
at various frequencies. The thick line is X'(0, T) adopted for 

fitting x(a',T). 

FIG. 5. Real and imaginary parts of the experimental di- 
electric susceptibility (thick lines) and fitting curves (thin 
lines). 
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FIG. 6. Recovery behavior of the NMR echo signal after a 
RF pulse sequence saturating the central (+1/2 <^ —1/2) 
transition for different temperatures (squares 600 K, dia- 
monds 179 K and stars 120 K). At high temperature the 
recovery is well described by Eq. (8) (solid line). In the 
low temperature range the recovery shows the change over to 
stretched exponentials. 



FIG. 7. Temperature dependence of the NMR relaxation 
rates 2W for different measuring frequencies. The lines rep- 
resent the theoretical curves assuming the LC distribution 
(see text). 



FIG. 8. Comparison between the NMR relaxation rates 
2W and the J{uj, T) curves extrapolated from the fits to 
x"(a>,T) of Fig. 5. 



FIG. 9. 2W relaxation rates at T = 140 K as a function of 
the measuring frequency (filled squares). The lines give the 
frequency dependence expected for the uniform distribution 
of the correlation time (dashed line) and for the LC model 
(dotted-dashed line). The solid line is obtained by adding to 
the LC behavior a frequency independent contribution of 0.29 



FIG. 10. Comparison between the imaginary parts of the 

dielectric susceptibilities (upper panel) and elastic compli- 
ances (lower panel) of PLZT 9/65/35 and PLZT 22/20/80. 
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